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ABSTRACT 

We present X-ray spectral analysis of the accreting young star TW Hydrae from a 489 ks observation using 
the Chandra High Energy Transmission Grating. The spectrum provides a rich set of diagnostics for electron 
temperature Tg, electron density A'^., hydrogen column density Nh, relative elemental abundances and velocities 
and reveals its source in 3 distinct regions of the stellar atmosphere: the stellar corona, the accretion shock, and 
a very large extended volume of warm postshock plasma. The presence of Mg XII, Si XIII, and Si XIV emission 
lines in the spectrum requires coronal structures at ^ 10 MK. Lower temperature lines (e.g., from O VIII, Ne IX, 
and Mg XI) formed at 2.5 MK appear more consistent with emission from an accretion shock. He-like Ne IX 
line ratio diagnostics indicate that = 2.50 ±0.25MK and A^^. = 3.0 ±0.2 x 10'^ cm"^ in the shock. These values 
agree well with standard magnetic accretion models. However, the Chandra observations significantly diverge 
from current model predictions for the postshock plasma. This gas is expected to cool radiatively, producing 
O VII as it flows into an increasingly dense stellar atmosphere. Surprisingly, O VII indicates Ne = 5.7^| 2 ^ 10" 
cm"^, five times lower than A^^. in the accretion shock itself, and ^ seven times lower than the model prediction. 
We estimate that the postshock region producing O VII has roughly 300 times larger volume, and 30 times more 
emitting mass than the shock itself. Apparently, the shocked plasma heats the surrounding stellar atmosphere 
to soft X-ray emitting temperatures and supplies this material to nearby large magnetic structures - which may 
be closed magnetic loops or open magnetic field leading to mass outflow. Our model explains the soft X-ray 
excess found in many accreting systems as well as the failure to observe high A,, signatures in some stars. Such 
accretion-fed coronae may be ubiquitous in the atmospheres of accreting young stars. 

Subject headings: accretion — stars: coronae — stars: formation — stars: individual (TW Hydrae) — tech- 
niques: spectroscopic — X-rays: stars 
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1. INTRODUCTION 

Low mass stars in star-forming regions produce strong X- 
ray emission from coronal magnetic activity, as evidenced by 
high temperature (^10 MK) emission from flares and active 
regions (Feigelson & Montmerle 1999; Gagne et al. 2004; 
Preibisch et al. 2005). The role of accretion in the production 
of X-rays from young stars is less well understood. Magnetic 
structures on stars can extend all the way to the inner accre- 
tion disk, as suggested from models of decaying stellar flares 
(Favata et al. 2005). Extended X-ray emission from jets also 
occurs, as observed in DG Tau, a system oriented so that the 
circumstellar disk plane aligns with our line of sight (Giidel 
et al. 2008). Soft X-ray emission may arise from shocks at 
the stellar base of this extended X-ray jet (Giidel et al. 2008; 
Schneider & Schmitt 2008). 

X-rays can also be produced in accreting systems when 
the accreting material accelerates to supersonic velocities and 
shocks near the stellar surface, heating the gas to a few MK. 
In the standard model of accreting Classical T Tauri Stars 
(CTTS), the magnetic field of the star connects to the accre- 
tion disk near the corotation radius, channeling the accretion 
stream from the disk to a small area on the star (Calvet & Gull- 
bring 1998). X-rays from the shock may be difficult to detect 
if the shock is formed too deep in the photosphere (Hartmann 
1998) or if they are absorbed by the stream of preshocked neu- 
tral or near-neutral gas. Furthermore, the X-ray shock signa- 
ture may be difficult to distinguish from soft coronal emission 
without the high resolution spectroscopy necessary to deter- 
mine temperature and density and identify the emitting re- 
gions. Shock models predict high electron density Ne (^ 10'^ 
cm"^) at relatively low electron temperature Te (a few MK), 
features that can be tested with X-ray fine ratio diagnostics. 



TW Hydrae (TW Hya) was the first CTTS found to show 
these characteristics of high density and low temperature 
(Kastner et al. 2002; hereafter, K02), based on a 48 ks Chan- 
dra High Energy Transmission Grating (HETG) spectrum. 
TW Hya is one of the oldest known stars (~' 10 Myr) still 
in the CTTS (accreting) phase. It is uniquely poised in an 
interesting state of evolution when it will soon stop accret- 
ing, lose its disk, and perhaps form planets (e.g., Calvet et al. 
2002). At a distance of only 57 pc, it is also one of the bright- 
est T Tauri stars in X-rays. Interstellar absorption is minimal, 
and the observed absorption can be considered intrinsic to the 
system (K02). TW Hya's circumstellar disk is nearly face-on 
with an inclination angle of 7° (Qi et al. 2004). Thus the star 
appears pole-on. 

Key to K02's argument for accretion is the electron density 
{Ne ^ 6 X 10'^ cm"-') determined from He-like Ne IX emission 
line ratios, and confirmed by subsequent observations with 
the XMM-Newton Reflection Grating Spectrometer (Stelzer 
& Schmitt 2004) and the Chandra Low Energy Transmission 
Grating (LETG; Raassen 2009). While A^^ for some active 
cool star coronae approaches such high values at high Te{^ 10 
MK), e.g., 44 Boo, UX Ari, and II Peg (Sanz-Forcada et al. 
2003; Testa et al. 2004), the high A^^ of TW Hya is produced at 
the significantly lower Te (~3 MK) expected for the accretion 
shock. BP Tau (Schmitt et al. 2005), V4046 Sag (Gunther et 
al. 2006), RU Lup (Robrade & Schmitt 2007), and MP Mus 
(Argiroffi et al. 2007) also show high values of the electron 
density. On the other hand, T Tau and AB Aur have a low 
Ne (Gudel et al. 2007a; Robrade & Schmitt 2007); and, the 
binary Hen 3-600 appears to lie at an intermediate value of A^^ 
(Huenemoerder et al. 2007). Despite large differences in A^^., 
all of these stars show an excess ratio of soft X-ray emission 
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to hard X-ray emission, manifested by their large O VII/O VIII 
ratios as compared with active main sequence stars (Giidel & 
Telleschi 2007; Robrade & Schmitt 2007). Additional soft X- 
ray diagnostics available with a deep Chandra exposure might 
be expected to shed light on the accretion process and test the 
shock models. 

Since accretion is expected to result in higher densities than 
found in an active stellar corona, it is critical that the density 
determination be secure. Interpretation of the Ne IX line ratio 
in terms of high density is not unique, since photoexcitation 
by ultraviolet radiation mimics the effects of high density. De- 
termining the density from Mg XI is not subject to the same 
degeneracy, since the relevant ultraviolet radiation, observed 
with the Far Ultraviolet Explorer (FUSE; Dupree et al. 2005), 
is too weak to photoexcite its metastable level. The three He- 
like systems (O VII, Ne IX, and Mg XI) together can provide 
independent determinations not only of the electron density, 
but also of the electron temperature and hydrogen column 
density. Deep spectroscopy can, first, put the accretion in- 
terpretation on firm ground and second, allow us to determine 
the structure of the accretion region itself. 

We were awarded a Chandra Large Observing Project to 
establish definitively whether the X-ray emission in TW Hya 
is associated with accretion, using high signal-to-noise ra- 
tio spectroscopy with the HETG, and, if so, to determine 
the emission measure distribution and elemental abundances, 
search for X-ray signatures of infall, outflow, and turbulence, 
and constrain the properties of the accretion hot spot. Dur- 
ing the Chandra observations we also conducted a ground- 
based observing campaign, which will be reported elsewhere 
(Dupree et al. 2010, in preparation). Here, we describe the 
X-ray observations in Section 2 and report analysis results 
in Section 3, where the densities, temperatures, and column 
densities can be extracted from this rich spectrum. Section 3 
also presents elemental abundances derived from fits to the 
spectrum, as well as the measurement of turbulent velocity. 
Section 4 describes spectral predictions from a model of the 
TW Hya shock. In Section 5 we discuss evidence for a new 
type of coronal structure. Section 6 gives a summary of results 
and conclusions. 

2. OBSERVATIONS AND SPECTRAL ANALYSIS 

Chandra observed TW Hya with the HETG in combination 
with the ACIS-S detector array for a total of 489.5 ks inter- 
mittently spanning the period between 2007 Feb 15 and 2007 
Mar 3. The observation consists of 4 segments, with expo- 
sure times as follows: 153.3 ks (Obsid 7435), 157.0 ks (Obsid 
7437), 158.4 ks (Obsid 7436), and 20.7 ks (Obsid 7438). The 
summed spectrum has about ten times the exposure time of 
the spectrum analyzed by K02. Figure 1 shows the total first- 
order light curve, with flux variations of about a factor two 
during the observation. Assuming a distance of 57 pc, the 
observed X-ray luminosity Lx, measured using the Medium 
Energy Grating between 2.0 and 27.5 A, is 1 .3 x 10^° erg s~' . 
One clear flare, with total X-ray luminosity of 2.1 x 10^° erg 
s"', persists for ^15 ks, producing a total energy of 4 x 10^"* 
ergs. The flare is retained in the following analysis since its 
contribution to line fluxes is negligible (see Section 3.4). 

We have extracted spectra from the event files, and pro- 
duced calibration files using the Chandra Interactive Analysis 
of Observations (CIAO ver 3.4) software package (Fruscione 
et al. 2006). With response matrices combined from the four 
observations, Sherpa (Freeman et al. 2001) was used to ex- 
tract individual emission line fluxes from the co-added obser- 
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Fig. 1 . — X-ray light curves (first order counts vs time since start) for the 
four observing segments of the Chandra observation. The hashed area repre- 
sents a gap of 650 ks. The upper curve shows the light curve for the total first 
order counts. This light curve includes both line and continuum emission. 
The middle curve uses the combined line counts from N VII, O VII, O VIII. 
Ne IX, Fe XVII, and Mg XI (identified as the "accretion" lines using Model 
D, discussed in Section 3.4). The lower curve uses the combined line counts 
from Mg XII, Si XIII, and Si XIV (identified as the lines from the "corona" 
from Model D). Note that the line-based light curves also include any contin- 
uum counts within the line profile. The data are binned over 3 ks. 

vations, maintaining separate first order spectra for the four 
grating arms (plus and minus; MEG and HEG) during the fit- 
ting procedure. A model for the continuum plus background 
was found from a global fit to line-free regions and used in ad- 
dition to the line models. Background rejection is high with 
ACIS order sorting. Thus the background contributes rela- 
tively little to this model except at the long wavelength end 
of the spectrum, where the effective area is low. In the O VII 
region, for example, the background contains about half as 
many counts as the continuum. For the relatively weak but 
important O VII line, background, continuum, and line have 
about 1, 2, and 16 counts, respectively. 

Line models are Gaussian functions with negligible or small 
widths relative to the width of the instrumental line response 
functions (-0.012 A and 0.023 A, for HEG and MEG, re- 
spectively). Table 1 gives the identified lines and their ob- 
served fluxes. For most lines the fit widths are less than half 
the instrumental width and are not statistically significant. We 
discuss the widths of the strongest lines in Section 3.5. 

All line fluxes are within 3<t of the values reported by K02 
except for Ne IX A13.437, which has slightly more than twice 
as much flux in our observation. Our line fluxes are also 
within 3ct of the values reported by Stelzer & Schmitt (2004) 
for lines observed with both instruments, except for Ne IX 
A13.553, for which our flux is about half. For lines in com- 
mon with those reported by Raassen (2009) all fluxes agree 
within 3cr. With the higher signal-to-noise ratio of our data, 
new diagnostic information is available. 

A major goal of this study is the determination of electron 
density from the He-like ions. We first discuss the diagnos- 
tics, and then model the global spectrum. Figure 2 shows the 
high quality in the He-like O VII and Ne IX line regions. In- 
terestingly, the Ne IX region shows no apparent Fe XIX lines, 
which can contaminate the Ne IX diagnostics (e.g., Ness et 
al. 2003). Upper limits on the strongest Fe XIX line at 13.52 
A imply that Fe XIX contributes 1 % or less to the Ne IX line 
fluxes. 

The Mg XI triplet line region is more complex (Fig. 3), 
with lower signal-to-noise ratio and apparent line blending 
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TABLE 1 
Selected Emission Lines 



Line 




h 

Flux 


Lit Area 


... 
Origin 




(A) 


(10 ^ ph cm ^ s ^) 


(cm"^) 




Si XIV 


6 182 


1 87 di 19 


163 




Si XIII 


6 648 


9 79 _i_ n n 

z. /z ^ VJ. i / 


146 6 


corons 


Si XIII 


6 688 


U. / 1 ^ W. it 


138 3 


coronfl 


Si XIII 


6 740 


1 .DO ^ W. i J 


1516 


coronfl 


MsT XTT 


8 421 


9 94 _i_ n 9n 


162 4 




Mq XI 


9.169 


2.25 ± 0.22 


135.4 


shock 


Mg XI 


9.231 


0.90 ± 0.17 


133.2 


shock 


Me XI 


9.314 


1.27 ± 0.19 


130.1 


shock 


NeX 


9.481 


1.64 ± 0.20 


114.5 


corona 


NeX 


9.708 


2.75 ± 0.27 


86.2 


corona 


NeX 


10.239 


9.23 ± 0.53 


81.1 


corona 


NelX^ 


10.765 


3.54 ± 0.43 


67.9 


shock 


Ne IX 


11.001 


8.48 ± 0.63 


67.4 


shock 


Ne IX 


11.544 


23.7 ± 0.99 


55.3 


shock 


Fe XXII 


11.783 


3.28 ± 0.50 


52.3 


corona 


Fe XXII 


11.932 


3.01 ± 0.50 


49.6 


corona 


NeX 


12.134 


76.1 ± 1.9 


48.5 


corona 


Fe XVII 


12.266 


4.3 ± 0.6 


46.6 


llllACU 


Ne IX 


13 447 




22 3 


shock 


Ne IX 


1 J.J J J 


114 ^ ^ 
i it. m J.J 


22 7 




Ne IX 


ij.oyy 


<o 7 ^ 9 'J 
Jo. / HI i. J 


25 


shock 


Fe XVIII 


14 208 


4 <Q -1- 7fi 

t. J^ u. / u 


29 




Fe XVIII 


14 256 


1 .\Jj W — J J 


28 5 




O VIII 


14 820 


s 04 -t- sn 


91 7 


shock 


Fe XVII 


15 014 




9 1 Q 




O VIII 


15 176 


10 4 i 1 2 


20 8 


shock 


Fe XVII 


15 261 


17.4 ±1.4 


20 


shock 


O VHP 


16.006 


30.6 ± 2.3 


15.9 


shock 


Fe XVIII 


16.071 


5.17 ± 1.0 


16.3 


corona 


Fe XVII 


16.780 


25.5 ± 3.0 


12.1 


shock 


Fe XVII 


17.051 


27.9 ± 2.3 


11.1 


shock 


Fe XVII 


17.096 


26.3 ± 2.3 


11.2 


shock 


O VII 


17.207 


4.68 ± 1.4 


10.7 


shock 


O VII 


17.396 


4.68 ± 1.5 


9.8 


shock 


O VII 


17.768 


3.66 ± 1.6 


7.7 


shock 


O VII 


18.627 


17.9 ± 1.4 


7.3 


shock 


O VIII 


18.969 


213. ± 8.4 


6.6 


shock 


N VII 


20.910 


3.8 ± 2.5 


3.3 


shock 


O VII 


21.601 


117. ± 10.0 


2.8 


postshock'' 


O VII 


21.804 


72.4 ±9.1 


2.7 


postshock'' 


O VII 


22.098 


15.2 ±4.4 


2.1 


postshock*" 


N VII 


24.781 


68.5 ±7.6 


2.6 


shock 



" Reference wavelength from ATOMDB. Wavelengths for multiplets are 
given as intensity-weighted averages. 

Observed fluxes at Earth, with no correction for absorption. Flux errors 
are la. 

Effective areas are estimate by combining the four grating arms and aver- 
aging the effective areas over the line width. 

The origin of the emission lines is based on Model D parameters (Section 
3.4). 

■= Blend with Fe XVII A10.77. 

' It is interesting to note that this Fe XVII line is mixed, while the longer 
wavelength lines are primarily from the shock. The reason for this differ- 
ence is that the 12.266 A line flux is relatively stronger at higher tempera- 
ture. 

s Blend with Fe XVIII A16.004; however, the predicted Fe XVIII contri- 
bution is negligible. 

This line is primarily formed in the postshock cooling plasma (Section 5). 
In Model D the EM distribution has only one component representing both 
the shock and postshock regions. 

between 9.2 and 9.4 A. For this region, the fit is constrained 
such that all known lines are fixed in wavelength relative to 
the Mg XI resonance line A9.169, which is fit separately for 
each grating arm. Lines are Gaussian functions with widths 
set to zero (i.e., only the instrumental line response function 
is used). Line fluxes of the Ne X Lyman series lines « = 5 
to 10 are scaled according to their oscillator strengths, since 
collision strengths are not available in the literature. We note 
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Fig. 2. — Top: Observed MEG spectrum in the O VII spectral region, 
shown as a histogram, with best-fit model (3 Gaussian lines) overlaid. Wave- 
length positions are fixed relative to each other and the instrumental line re- 
sponse function is included in the model. Resonance (r), intercombination 
(i) and forbidden (f) lines are marked. Bottom: Observed MEG (upper) and 
HEG (lower) spectra in the Ne IX region, shown as histograms, with best- 
fit model overlaid. For this fit, the positions, fluxes and line widths are free 
parameters, yielding the constraints on velocities described in the text. 



that systematic errors from line blends are not included in the 
errors given in Table 1 . 

3. RESULTS 

One of the advantages of X-ray spectroscopy with high 
signal-to-noise characteristics is the availability of numer- 
ous line ratio diagnostics to assess conditions in the radiat- 
ing plasma. In particular. He-like systems provide valuable 
diagnostic opportunity, as discussed extensively in the liter- 
ature (Porquet et al. 2001; Smith et al. 2001; and references 
therein). The G-ratio, defined as the flux ratio of the forbid- 
den {ls2s ^Si — > li^ '5o) plus intercombination (li2/? — > 
Is^ '5o) lines to the resonance (ls2p Is^ ^Sq) line, is 

sensitive to the electron temperature (Tg), while the R-ratio 
(flux ratio of the forbidden to intercombination lines) is sen- 
sitive to Ne- Ratios of emission lines in the resonance line 
series (transitions from 2p, 3p, etc. to ground, also known as 
Hcq', He^, and so on) are sensitive both to T^, through their 
relative Boltzmann factors, and to the absorption by an inter- 
vening hydrogen column density {Nh) through their energy 
separation. 



Hcrt is the resonance line sometimes labeled r or w. 
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9.15 9.20 9.25 9.30 9.35 
Wavelength (Angstroms) 

Fig. 3. — Observed MEG (upper) and HEG (lower) spectra in the Mg XI 
spectral region, shown as histograms, with best-fit model spectrum, described 
in text, overlaid. Positions of resonance (r), intercombination (i), and forbid- 
den (f) lines are marked. 

This analysis emphasizes the use of line ratio diagnostics, 
including those discussed above. We use the Astrophysi- 
cal Plasma Emission Code (APEC) with the atomic data in 
ATOMDB vl.32 (Smith et al. 2001) to calculate Hne emis- 
sivities, except for He-like Ne IX, for which we use the more 
accurate calculations of Chen et al. (2006). The model for ab- 
sorption by gas with cosmic abundances is taken from Morri- 
son & McCammon (1983). 

In this section we present the physical conditions of the 
emitting plasma as determined from the line ratios of specific 
He-like and other ions, in particular Ne, Tg, and Nh- We use 
information determined from line ratios to constrain an em- 
pirical model of the emission measure distribution (EMD), 
and then allow additional parameters of the model to be con- 
strained by global fits to the spectrum. We give the elemental 
abundances determined from these fits. Finally we present 
velocity constraints from line centroids and widths. 

3.1. Electron temperature 

We use G-ratios from multiple ions to determine the elec- 
tron temperature. Figure 4 shows the observed G-ratios of He- 
like O, Ne and Mg, overplotted on their respective theoretical 
functions of T^. For comparison, the emissivities of their reso- 
nance lines peak at Tg = 2, 4, and 6 MK, respectively, as shown 
in the figure. Table 2 gives the derived Tg from these ions. The 
Ne IX G-ratio indicates that = 2.5 ± 0.25 MK. The Ne IX 7; 
determination is highly reliable, given that recent theoretical 
calculations, which include the resonance contributions to the 
collisional excitation cross sections, agree with experimental 
measurements to within 7% (Chen et al. 2006; Smith et al. 
2009a, 2009b). On the other hand, the theoretical calculations 
for O VII and Mg XI have higher systematic uncertainties, and 
the observational errors are larger, such that the values de- 
rived from O VII and Mg XI carry less weight (but are consis- 
tent with Ne IX). Both Ne IX and Mg XI are found far below 
their temperature of maximum emissivity T„,ax, supporting an 
accretion shock scenario for their formation. 

3.2. Electron density 

We use R-ratios from multiple ions to determine the elec- 
tron density. Figure 5 shows the observed R-ratios for O VII, 
Ne IX, and Mg XI, overplotted on the theoretical functions of 

^ http://cxc.harvai'd.edu/atomdb 
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Fig. 4. — Theoretical G-ratios (i.e., ratios of forbidden plus intercombina- 
tion to resonance line fluxes) as functions of Te, for O VII (dashed), Ne IX 
(solid), and Mg XI (dash-dotted). Overplotted are the observed ratios with 
l(j errors. The model for Ne IX is from Chen et al. (2006). The models for 
O VII and Mg XI are from APEC using ATOMDB vl.3 (Smith et al. 2001). 
Scaled emissivities (ph cm^ s"') for the resonance lines of the three ions are 
also shown. For Ne IX and Mg XI the derived temperatures are well below 
the temperatures of their peak emissivities, consistent with an intepretation 
of their formation in an accretion shock. 



TABLE 2 
He-like G-Ratio Te Diagnostics 



Ion 


Observed Ratio" 


(K) 


OVII 
NelX 
MgXI 


0.75±0.13 
0.98±0.04 
0.97±0.19 


2.4!^-^ X 10*^ 
2.5!<]|5 X lO'S 
2.1+2-; X 10*5 



' Errors are Icr. 



TABLE 3 
He-like R-RatioW, Diagnostics 



Ion 


Observed Ratio" 


Ne 






(cm-^) 


OVII 


0.21 ±0.07 


5.7![^ X 10" 


NelX 


0.51 ±0.03 


3.0+g| X 10'2 


MgXI 


1.41±0.34 


X 10'2 



" Errors are Icr. 

electron density A^^, from APEC. Table 3 presents the Ne val- 
ues derived from these R-ratios. Sensitivity to Te is negligible 
near the observed R-ratio values. Smith et al. (2009) show 
that the APEC R-ratio for Ne IX in the high density range is 
indistinguishable from the R-ratio calculated by Chen et al. 
(2006). The theoretical R-ratios from APEC for all three ions 
are in excellent agreement with Porquet et al. (2001) over the 
density range reported here. 

From Ne IX we derive A^^ = 3.0 ±0.2 x 10'^ cm"^ while 
O VII gives Ne = 5.1^^2 x 10" cm'l The 3a range from 
Ne IX, between 2.5 and 3.9 x 10'^ cm"-*, provides a tight con- 
straint.^ The best-fit value for Mg XI, using all the lines in the 
region as described in Section 2, is A^^ = 5.8;I;24 x 10'^ cm""*. 

^ The measured value of the Ne IX R-ratio is 0.515 ± 0.033, quite close to 
the value of 0.446 ± 0. 124 reported by K02. The difference of nearly a factor 
of 2 in the derived Ne appears to be due to the difference in atomic data used, 
with their models possibly based on the Raymond & Smith code (Raymond 
1988). 
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Fig. 5. — Theoretical R-ratios (i.e., ratios of forbidden to intercombination 
line fluxes) as functions of Ne for He-like O, Ne and Mg. Overplotted are 
the observed ratios with Icr errors. Models are from APEC using ATOMDB 
vl .3. While Figure 4 shows that these ions are formed at similar temperature, 
they are found at significantly different densities. 



For comparison, we have also fit only the 3 lines of Mg XI, 
again with line widths and relative positions fixed, obtaining 
Ne = l ±2 X 10'^ cm"-', in good agreement with the best-fit 
value. This comparison suggests that line blending does not 
unduly affect the result. 

At low Ne the population of the metastable 2s level 
builds up since the decay rate is extremely low. At high A^^ 
electron impact excites the metastable level population up to 
the 2p^P levels. In principle, low R-ratios may also be pro- 
duced at low Ne by photoexcitation. Stelzer & Schmitt (2004) 
and K02 dismiss photoexcitation based on the stellar temper- 
ature 4000 K), but Drake (2005) suggests that a hot spot 
above the stellar surface, where the accretion stream shocks, 
could produce sufficient photoexciting radiation. We have 
measured the flux at the excitation wavelengths 1647 A, 1270 
A and 1034 A for O VII, Ne IX, and Mg XI, respectively, from 
the Hubble Space Telescope STIS (Herczeg et al. 2002) and 
the FUSE specti-a (Dupree et al. 2005). The intensity at 1034 
A is an order of magnitude lower than at the longer wave- 
lengths, constraining the black body temperature of the hot 
spot to be less than 10,000 K, in good agreement with Inter- 
national Ultraviolet Explorer analysis showing a blackbody 
at 7900 K covering about 5% of the stellar surface (Costa et 
al. 2000). Thus the ultraviolet measurements and the Mg XI 
R-ratio rule out photoexcitation and support the interpretation 
of high A^,,. Together with the temperature derived from the 
Ne IX G-ratio, the density supports an accretion origin for 
Ne IX. 

The ratio of Fe XVII A 17. 096 to A 17. 051 becomes sensitive 
to Ne above 10'^ cm"' (Mauche et al. 2001). The mea- 
sured ratio of 1 .06 ± 0. 1 3 is at the low A^^ limit, ruling out the 
high density from this line ratio suggested by Ness & Schmitt 
(2005). 

3.3. Hydrogen column density 

The HETG spectrum contains several series of resonance 
lines that indicate absorption by neutral gas along the line 
of sight. The most useful of these diagnostics are the O VII 
and Ne IX He^/He^ ratios because their formation tempera- 
tures are independently determined from the G-ratio. Table 4 
presents the observed line ratios, and their use as diagnos- 
tics for Nh- Figure 6 shows the theoretical He^/He^ ratios 
as functions of Te for O VII and Ne IX, for both unabsorbed 



TABLE 4 
Line Ratio A'h Diagnostics 



Identification 


Observed 


Predicted Ratio for log 


Nh (cm-2): 




Ratio" 


0.00 


20.61 


20.89 


21.26 


O VII Heo/He/?'' 


6.54±0.76 


8.54 


6.52 


5.22 


2.49 


VIII Lya/Ly/3'' 


6.95±0.58 


9.63 


8.07 


6.82 


4.55 


Ne IX Hea/He/3'' 


7.38±0.35 


11.14 


10.07 


9.23 


7.34 


Ne X Ly«/Ly^'' 


8.24±0.52 


12.93 


11.89 


10.91 


8.93 


Ne X Lya/Ly/S"^ 


8.24±0.52 


7.061 


6.49 


5.96 


4.88 



" EiTors are Ict. 
Predictions assume Te = 2.5 X lO'' K. 
Predictions assume = l.Ox 10^ K. 




o r 1 1 1 1 1 

6.1 6.2 6.3 6.4 6.5 

log [Te (K)] 

Fig. 6. — Theoretical models of Hcq to Heg line flux ratios as functions of 
Tf, with solid curve for Ne IX and dashed curve for O VII. The filled circle 
is the observed ratio for Ne IX and the open circle is for O VII, with Icr 
error bars. Unabsorbed and best absorbed models for each ion are labeled. 
The absorbed models shown have A'a = 1.8 X 10"' and 4.1 X 10^'' cm"^ for 
Ne IX and O VII, respectively. 

and absorbed cases. The observed ratios are overplotted us- 
ing Te determined from their G-ratios (which, as noted above, 
is more reliable for Ne IX than for O VII). For the best-fit 
Te, Ne IX gives Nh = 1.8 ±0.2 x 10^' cm"^ and O VII gives 
Nh = 4.l!!![ g X 10^" cm"^. These two values show a mean- 
ingful difference. For Ne IX at 2.5 MK, the observed ratio is 
1 1(7 below the unabsorbed model, making the presence of an 
absorbing column quite definitive. The column density from 
O VII is consistent with the lower limit of Nh = 2 x 10^° cm"^ 
determined from the absence of Ne VIII and Si XII lines long- 
ward of 40 A in the Chandra LETG spectrum (Raassen 2009). 

Ratios of the H-like series lines (Ly^, Ly^, etc.) are also 
sensitive to both absorption and temperature, but unlike the 
He-like ions, H-like ions do not present an independent tem- 
perature diagnostic. If the absorption column depends on the 
ion, as implied by the inconsistency between Ne IX and O VII, 
then we need to assume a model Tg to obtain a column density 
from the H-like line ratios. Assuming 2.5 MK, the observed 
ratio O VIII Ly^/Ly^ indicates Nh = 7.8!^;^ x 10^" cm"^. This 
value is consistent within the errors with Nh from O VII and 
thus a third absorber is not required. For Ne X at 2.5 MK, 
the column density would be larger than the column density 
for Ne IX, as shown in Table 4; however. Table 4 also shows 
that no absorption is required if the Ne X lines are formed at 
10 MK. We must further consider the effects of temperature 
on Ne X line ratios. 

We can also estimate the temperature using the charge state 
balance. Table 5 presents the Hea/Lyo, ratios for oxygen and 
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TABLE 5 
Te From Charge State Balance 



Identification 


Observed^ 


Corrected for log 


Nh (cm-2): 






20.61 20.89 


21.26 


O VII Hen/O VIII Lvn 


0.55±0.052 


0.70 0.85 


1.67 


Nc IX Hca/iXc X Lya 


2.32±0.()8 


2.48 2.60 


3.04 


Identification 




Derived (10* K)'' 




OVIlHea/0 VniLya 


3.03 


2.83 2.66 


2.23 


Ne IX Hea/Ne X Lya 


4.16 


4.08 4.01 


3.87 



" Errors are \a. 

The derived Te values are based on the ratios listed above in the corre- 
sponding observed and corrected columns. 



neon. First, the observed ratios are corrected for different col- 
umn densities as found above, and then Te is derived from 
the respective corrected ratios. For oxygen, good consistency 
is found for the temperatures derived from the G-ratio and 
the O vn Hea/O VIII Ly^ ratio, given a column density of 
about 8 xlO^" cm"^. For neon, on the other hand, all the de- 
rived temperatures are larger than T^. derived above from the 
Ne IX G-ratio. Thus it appears that some of the Ne X emis- 
sion originates at higher temperature. Using Nh = 2 x 10^^ 
cm~^ and Tg = 2.5 MK, and the observed Ne IX He^ flux, we 
predict the flux of the Ne X Ly^ to be only 18% of the ob- 
served value. In this case, the rest of the Ne X emission must 
be formed at higher temperature. If the trend of increasing 
Nh with charge state continues, then Nh for Ne X would be 
higher than 2 x lO'^' cm"^, and the fraction of Ne X formed at 
2.5 MK would be somewhat larger. 

The emissivity functions for lines of Fe XVII and Ne X 
both peak near 5 MK, suggesting the possibility of comparing 
their formation temperatures. The ratio of Fe XVII A15.014 
to A 15. 261 is sensitive to Tg, due to a blend of an iimer sheU 
excitation line of Fe XVI at 15.261 A (Brown et al. 2001). 
The temperature sensitivity essentially comes from the charge 
state balance of Fe XVI and Fe XVII. Using the reported ex- 
perimental dependence of the line ratio on the electron beam 
energy of the Electron Beam Ion Trap at Lawrence Livermore 
National Laboratory (Brown et al. 2001), the observed HETG 
value gives = 2.34;I;y J]^ MK. Despite similar peak temper- 
atures, Fe XVII and Ne X are not formed together. Appar- 
ently, the long tail of the Ne X Ly^ emissivity function to- 
ward high temperature contributes significantly to its emis- 
sion. The Fe XVII lines are predominantly formed in the ac- 
cretion shock, while the Ne X lines are formed both in the 
shock and in the hotter corona. 

One might expect that the ions formed at low temperature 
(O VII, O VIII, and Ne IX) should all experience similar ab- 
sorption. Using Nh derived from Ne IX, we can derive char- 
acteristic emitting Tg for O VIII and O VII of 1.4 MK and 
< 1.0 MK, respectively. These temperatures are far lower 
than Te derived from the O VII G-ratios. They are both also 
far lower than Te>2 MK derived from the charge state ratio 
of oxygen represented over the range of Nh found here, as 
shown in Table 5. We conclude that an increase in Nh with 
increasing charge state (presumably temperature) up through 
Ne IX is consistent with all the data. 

We note that Argiroffi et al. (2009) claim evidence for reso- 
nance scattering of the strongest lines in the XMM-RGS spec- 
trum of MP Mus, though their RGS spectrum of TW Hya is 
too weak to show the signatures of absorption we find here. 
We rule out resonance scattering in TW Hya using the strong 



series of Ne IX. An optical depth r of about 0.4 for Ne IX 
Hca would be required to model the observed Hca/He^ ra- 
tio. For a single ion, the optical depth t scales as gfoscK such 
that we can then predict the other Ne IX line ratios. We find 
then that Hca/He-y ratio is overpredicted by 50%. Since the 
oscillator strengths decrease rapidly with increasing princi- 
pal quantum number n, only the He^ line might be affected 
by resonance scattering. We find instead, that all the Ne IX 
Unes are affected by absorption, and rule out resonance scat- 
tering as the absorption mechanism. We conclude that the ab- 
sorber is neutral or near neutral, consistent with the preshock 
accreting gas, as suggested by theoretical studies. For exam- 
ple, Lamzin (1999) considers the absorption by the preshock 
gas of the X-ray emitting plasma for different geometries and 
orientations. Gregory et al. (2007) use reaUstic coronal mag- 
netic fields coupled with a radiative transfer code to calculate 
the obscuration of X-ray emission by accretion columns, and 
suggest that this effect can explain the observed low X-ray lu- 
minosities of accreting young stars relative to non-accretors. 

3.4. Emission measure (EM) distribution and elemental 
abundances 

We present here four models for the emission measure 
(EM) distribution"^ and elemental abundances (Table 6). 
These models serve to illustrate how values derived from the 
line ratios affect the global fit to the spectrum and to show 
whether abundance determinations are robust. For all four 
models, the first-order HEG and MEG spectra are fit to a set 
of variable abundance APEC models, using Sherpa and apply- 
ing different constraints to the fitted parameters. All models 
have acceptable goodness-of-fit values. Model A has two Te 
components and a single absorber. Model A should be appro- 
priate for comparison with other X-ray spectra of cool stars 
obtained at lower spectral resolution. 

Model B is also a two- 7], model, but with constraints im- 
posed from the line ratios. The lower T^ is fixed at 2.5 MK, 
and the single absorber Nh is fixed at 1.0 xlO^^ cm~^ (a 
rough average of the values obtained from line ratios). The 
A'e- sensitive forbidden and intercombination line regions are 
excluded from the fit. 

Table 6 also gives results from multicomponent Models C 
and D, each with the low Te and a single A^^ fixed, as in 
Model B. Assuming that the hotter component is coronal in 
nature, a broad distribution of Tg is reasonable, although the 
shape is not well constrained. We broaden the hot component 
using a Gaussian-shaped function centered around 12.5 MK. 
For Model C the entire spectrum is fit simultaneously, with 
only the forbidden and intercombination line regions excluded 
as for Model B. 

For Model D the emission measures for the two compo- 
nents are fit to line-free regions, identified visually and using 
APEC models. The thermal continuum emission fit to line- 
free regions is strongly temperature-dependent and thus con- 
strains the emission measure distribution. The abundances are 
then fit to narrow regions centered on the strong Unes. These 
separate fitting procedures are iterated until the values stop 
changing. The resulting Model D EM distribution is shown 

^ Emission measure EM is defined as J O.^N^dV , where the factor 0.8 
accounts for the hydgrogen to electron density ratio and the integral is taken 
over the volume V. The intensity of a spectral feature I = e EM /(AttD^) 
where e is the emissivity in units of ph cm^ s"' and D is the distance to the 
source. The function s depends on Te, and is given in steps of log [Te (K)] = 
0.1 in ATOMDB (Smith et al. 2001). 
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TABLE 6 

Model Parameters" from Global Spectral Fitting 



Parameter 


Model A 


Model B 


Model C 


Model D 


ri (MK) 


3.58±0.04 


2.5* 


2.5* 


2.5* 


EMI (10" cm-3) 


1.11±0.13 


1.49±0.09 


1.56±0.10 


4.22±0.70 


Tl (MK) 


18.8±0.5 


11.2±0.2 


12.6^^ 


12.6'' 


EM2 (10^3 cm-3) 


0.33±0.01 


0.47±0.01 


0.121±0.004'= 


0.101±0.005<' 


Nh (10" cm-2) 


0.15±0.22 


l.O* 


1.0* 


1.0* 


N (8.05)^ 


0.66±0.14 


0.28±0.06 


0.53±0.10 


0.20±0.04 


(8.93)^ 


0.23±0.02 


0.14±0.01 


0.24±0.02 


0.09±0.01 


Ne (8.09f 
Mg (7.58)« 
Si (7.55f 
Fe (7.67)^ 


1.23±0.08 


2.04±0.11 


2.50±0.13 


1.04±0.02 


0.18±0.02 


0.16±0.02 


0.21 ±0.02 


0.18±0.02 


0.33±0.03 


0.24±0.02 


0.29±0.03 


0.34±0.03 


0.13±0.01 


O.lOiO.Ol 


0.20±0.01 


0.16±0.01 



" AH errors given are statistical errors from the fit. 
* Value is fixed. 

EM2 given is the fit to the peak of the EM distribution. Shape of the hot compo- 
nent is the same as for Model D, as described in the text. 

^ EM2 given is the fit to the peak of the EM distribution. Shape of the hot compo- 
nent EM distribution is shown in Figure 7 and described in the text. 

Abundance of element relative to the solar values of Anders & Grevesse (1989). 
The abundances in logarithmic form with H=12.0 are given in parenthesis next to 

the element label 




6.4 6.6 6.8 7.0 

log [T3 (K)] 

Fig. 7. — EM distribution for Model D as described in text. 



in Figure 7. The larger statistical error for the low T^. emis- 
sion measure in Model D compared with Model C (Table 6) 
stems from the smaller number of bins used in the fits; how- 
ever, we expect the systematic errors for Model D to be lower, 
given that we are only using bins whose information content 
is secure, and thus we choose Model D to illustrate the model 
spectrum.^ Figure 8 compares the Model D spectra predicted 
by the broad "coronal" component and the soft "accretion" 
component with the observed spectra. Table 1 lists the origin 
of each line as accretion or corona based on Model D. An ad- 
ditional component with Te lower than 2.5 MK can be added 
to both Models C and D but is not required, as it does not 
improve the fit significantly. Similar results are obtained with 
somewhat different choices of width and peak Tg for the hot 
component. 

Since the line ratios indicate increasing Nh with higher 
charge state, we tried adding a second absorber to the high 

' The four models are all statistically acceptable; however, it is interesting 
to note that none of the models do a good job of fitting all the emission line 
fluxes. For example, Model D does underpredicts Ne X Ijja by almost a 
factor of 2. We attribute these difficulties primarily to the complexity of the 
absorption, and secondarily to the few constraints on the shape of the EM 
distribution. 
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Observation 



lujM-Ajljlsi 



5 10 15 20 

Wavelength (Angstroms) 

Fig. 8. — Comparison of predicted and observed spectrum from the 
Medium Energy Grating (MEG). Upper curve is the predicted spectrum from 
the Model D EM distribution above 5 MK (the "coronal model"), shown in 
Figure 7. Middle curve is the predicted spectrum from the Model D emission 
measure at 2.5 MK (the "shock model"). The lower curve is the observed 
spectrum. As noted in the text. Model D significantly underpredicts the Ne X 
Lye . Also note that the forbidden and intercombination line predictions using 
the APEC models in Sherpa are at their low density values (e.g., Ne IX). 



Te Gaussian component for Models B, C, and D. Robust so- 
lutions are difficult to obtain, as the absorption, abundances, 
and normalizations of the emission measures do not indepen- 
dently determine the global X-ray spectrum. We also explored 
the possibility that the low and high Te components have dif- 
ferent metal abundances but again, without robust results. 

Using Te=2.5 MK, as derived from the Ne IX line ratio, in- 
stead of re=3.58 MK as found in Model A, has important im- 
plications for interpreting the formation of some of the emis- 
sion lines. The ratio of the strong Ne IX Hca and Ne X Ly„ 
line emission drives the Model A fit to more than 3.0 MK with 
very little absorption. Specifically, only ^15% of the Ne X 
Lya emission arises from the 2.5 MK component, whereas a 
3.58 MK component can produce more than half of the Ne X 
emission. The 2.5 MK component produces essentially all of 
the N VII, O VII, O VIII, and Ne IX emission and more than 
half the emission from Fe XVII. The fluxes from these lines 
formed at 2.5 MK show no increase over their average value 
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during the flare, as shown in Figure 1, indicating that the flare 
is associated with the hotter corona rather than the accretion 
shock, and justifying its retention in our analysis. Lines from 
Mg xn. Si Xin, Si xrv (Fig. l), and Fe XXII are produced by 
the hotter component and may participate in the flare. While 
K02 found that Mg XI was overpredicted by their model by a 
factor of about 3, our Models B, C, and D have no such prob- 
lem. In fact, for Models C and D, the small flux of the Mg XI 
resonance line forces the emission measure to be negligible 
between ~3 and 5 MK (Fig. 7). 

The relative elemental abundances given in Table 6 are rea- 
sonably consistent for the different models. The abundances 
are also similar to those found by K02, Stelzer & Schmitt 
(2004), and Drake et al. (2005). It is important to note that 
the continuum emission arises primarily from the hot compo- 
nent, except at the longest wavelengths (Fig. 8). Thus ab- 
solute abundances in the low component for Ne, O, Fe 
and other species cannot be determined reliably. Neverthe- 
less, the extremely large Ne/O abundance ratio appears to be 
a very robust result. Differences in the absolute abundances 
derived here reflect both differences in Tg and the degeneracy 
between emission measure and Nh- For example, nitrogen 
and oxygen are formed entirely in the low temperature com- 
ponent and thus the abundance differences among the models 
reflect differences among the low temperature emission mea- 
sure and Nh, as does the difference in the neon abundance 
between Models A and B. The difference in the neon abun- 
dance between Models C and D also reflects the differences 
in their emission measures at 2.5 MK. 

Weak emission is apparent from the S and Ar complexes, 
but we have not included line fluxes from these elements in 
Table 1 because the individual lines (in particular the He- 
like resonance line) cannot not be cleanly determined in the 
four grating arms. Instead we have used the Model C global 
fit to obtain 0.21 ± 0.09 and 1.41 ± 0.49 times solar (An- 
ders & Grevesse 1989) for the S and Ar abundances, respec- 
tively. A high Ar/O abundance would add weight to Drake et 
al.'s (2005) suggestion that the accretion stream is depleted of 
grain-forming elements; however, the high Ar abundance is 
unfortunately not statistically significant. 

3.5. Velocity constraints 

In this section, we discuss velocity constraints from the line 
measurements. The measurements of the emission line fluxes 
given in Table 1 use Gaussian lines in addition to the stan- 
dard HETG calibration line spread function. Table 7 gives the 
lines with the highest signal-to-noise ratios, for which Gaus- 
sian widths are measured. Lines not listed do not show sta- 
tistically significant widths, but are consistent with the range 
of widths found here. For these fits the continuum level is al- 
lowed to vary to provide a better local fit. Assuming the Ne IX 
lines form in the same region, and using a flat continuum flux 
level from 13.45 to 13.75 A, we force the 3 triplet lines to 
have the same Gaussian width to improve the significance of 
the velocity measurement. Figure 2 shows the best fit for the 
Ne IX lines. As before, the centroids of the lines measured 
in each grating arm are fit independently. We thus obtain a 
velocity of 183 ± 16 km s"'. The thermal Doppler FWHM 
velocity for neon at 2.5 MK is 75 km s"'. The error in the 
caUbration is estimated to be about 5%. Thus we determine 
a turbulent velocity of 165 ± 18 km s~\ including calibra- 
tion error in quadrature with the la statistical error. The line 
profiles are consistent with each other within the errors and 
suggest a range of velocities from zero up to about the sound 



speed of the gas (185 km s~' at 2.5 MK). The turbulent ve- 
locity is well below the preshock gas velocity of ^500 km 
s"^ Line centroids are consistent with reference wavelengths 
from ATOMDB to within the calibration uncertainty on the 
absolute wavelength scale. 

3.6. Physical conditions of the X-ray plasma 

We summarize here the physical conditions of the X-ray 
plasma from a consistent analysis of the spectra. The new di- 
agnostic measurements strongly support K02's argument that 
the low temperature X-ray component (Je ^ 2.5 MK) arises in 
the accretion shock. While the density we derive from Ne IX, 
Ne = 3 X 10^^ cm"^, is within the range found in active stars 
at higher Tg (^ 6 MK), it is at least an order of magnitude 
larger than the few Ne IX R-ratios reported (Huenemoerder et 
al. 2001; Ness et al. 2002; Osten et al. 2003). The A^^ derived 
from the O VII R-ratio is also more than an order of magni- 
tude larger than derived from O VII in other cool stars (see 
Sanz-Forcada et al. 2003; Ness et al. 2004; Testa et al. 2004). 
We also measure the Mg XI R-ratio and rule out photoexcita- 
tion, since the observed Mg XI R-ratio can only be produced 
by high A^,., namely 5.8 x lO'^ cm""*. 

The G-ratios from O VII, Ne IX, and Mg XI give similar 
Te ^2.5 MK, significantly lower than the Te of peak emis- 
sivity for Ne IX (4 MK) and Mg XI (6.3 MK). The G-ratio 
from Ne IX is particularly reliable because of good statistics 
and accurate atomic data. In our multi-Te models (Models C 
and D), this low component is isolated and does not con- 
nect continuously with the hotter, coronal component. Active 
main-sequence stars tend to have emission measures rising up 
to 6 MK and above, such that, if biased at all, Tg-sensitive line 
ratios should be biased toward higher rather than lower tem- 
perature. Thus the measurements also support the accre- 
tion shock scenario. More accurate G-ratio models for O VII 
and Mg XI are required to establish any differences among 
these ions. 

We also measure strong absorption using line ratios from 
several ions, and require at least two different absorbing col- 
umn densities, with the higher charge states experiencing 
more absorption. It seems Ukely that the gap between the two 
components of the EM distribution is caused by neutral H ab- 
sorption of the soft X-ray spectrum that is produced by the hot 
coronal component. The pattern of absorption is reminiscent 
of the "two-absorber X-ray sources" reported by Gtidel et al. 
(2007b), where the coronal component has ten times more ab- 
sorption than the soft component; however, unlike TW Hya, 
these sources have high accretion rates and are believed to 
drive jets. We note that the absorption in TW Hya is identi- 
fied here from line ratio diagnostics, whereas the absorption in 
Giidel et al. (2007b) was determined from CCD spectra. We 
were not able to determine multiple absorbers using standard 
global fitting methods, even with the high resolution spec- 
trum presented here, because the continuum spectrum has low 
signal-to-noise ratio. On the other hand, it seems likely that 
with low spectral resolution, multiple absorbers could also be 
easily missed. Thus the presence of more than one absorber 
may be a more universal characteristic than we can presently 
establish without reliable line ratio diagnostics. 

All the values of Nh found here are higher than those found 
using the hydrogen Lya profile, for which a conservative up- 
per limit is 6 X 10^^ cm~^ (Herczeg et al. 2004), indicating that 
the absorption is not due to the interstellar medium but is in- 
trinsic to the stellar system. Our column densities are consis- 
tent with the range of values previously reported from ROSAT 
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TABLE 7 

Velocities from Strong Emission Lines 



Line 


Kef 

(A) 


FWHM 


(km s~' ) 


(km s~' ) 


(km s~' ) 


NeX 


10.239 


0.0042 ± 0.0030 


122 ± 87 


75'' 


96 


NelX 


11.544 


0.0024 ± 0.0024 


62 ±62 


75 





NeX 


12.134 


0.0082 ± 0.0009 


203 ± 23 


75" 


189 


NelX 


13.447 


0.0103 ± 0.0010 


229 ± 22 


75 


216 


NelX 


13.553 


0.0057 ± 0.0014 


126 ± 32 


75 


101 


NelX 


13.699 


0.0073 ± 0.0019 


159 ±41 


75 


140 


Fe XVII 


15.014 


0.0065 ± 0.0040 


129 ± 80 


45 


121 


VIII 


18.969 


0.0074 ± 0.0026 


116 ±41 


85 


79 



^ Voi,j is the observed FWHM velocity. 

* Y,h is the thermal FWHM velocity calculated at 2.5 MK. 



Vf„rfc is the turbulent FWHM velocity assuming V,„rfc = y 
Note the thermal velocity for neon at 10 MK is 150 km s"' . 



and ASCA (Kastner et al. 1999). In fact the two different val- 
ues from ROSAT and ASCA are consistent with our different 
values and may reflect the complex absorption coupled with 
the different instrument responses rather than time-variable 
absorption. As noted also by Giidel et al. (2007b) for the 
highly absorbed coronae of DG Tau A, GV Tau, DP Tau, our 
X-ray measurements are higher than expected from the op- 
tical extinction assuming standard gas-to-dust ratios. For X- 
ray emission that is highly localized compared with the stellar 
Lya, and preferentially absorbed by the accretion stream di- 
rectly above it, a difference in Nh derived by the two methods 
seems entirely reasonable. 

4. THE ACCRETION SHOCK MODEL 

We can test the hypothesis that some of the measured X-ray 
emission comes from the accretion shock by constructing a 
standard one-dimensional (ID) model of the magnetospheric 
infall, shock heating, and postshock cooUng. This model is a 
simplified version of similar ID models in the Uterature (Cal- 
vet & Gullbring 1998; Gunther et al. 2007), and is tailored to 
specific measured properties of TW Hya. 

The magnetospheric accretion is assumed to follow a set 
of dipole magnetic field fines that thread the accretion disk. 
We use Konigl's (1991) expression for the inner truncation 
radius of the disk to determine the distance at which parcels 
of gas are launched from rest. This expression depends on the 
mass accretion rate Mace, the surface magnetic field strength 
B*, and the stellar mass M* and radius For TW Hya, 
we take M« = 0.7Mq and /?* = O.SRq (Batalha et al. 2002). 
As in Cranmer (2008), we assume a canonical T Tauri star 
magnetic field strength of « 1000 G. A range of accretion 
rates for TW Hya has been reported, from 4 x 1O~^''M0 yr"' 
(Muzerolle et al. 2000) to 5 x IQ-'^ Mq yr"' (Batalha et al. 
2002). The low end of this range produces the best agreement 
with the observations presented here, so the Muzerolle et al. 
accretion rate will be adopted for the cooling-zone models 
below. 

With the above values, the Konigl (1991) truncation radius 
of the disk is found to be approximately = 4.5/?*, which 
impHes a ballistic free-fall velocity at the stellar surface of 

-,1/2 



Vff = 



2GM* 



^ R* 
r, 



■ 509 km s" 



(1) 



Since the preshock gas is often assumed to have a sound speed 
only of order ~ 10 km s~' (e.g., Muzerolle et al. 2001), the 



highly supersonic flow should produce a strong shock at the 
stellar surface, heating the plasma to a postshock temperature 
of Ijjost ~ 3mv^/16k, where k is the Boltzmann constant and m 
is the mean atomic mass. For the above parameters, this value 
is Tposi = 3.4 MK. [We neglect any preshock heating due to 
photoionization since the X-ray luminosity is relatively low. 
See Calvet & Gullbring (1998).] 

The preshock mass density ppre in the accretion stream can 
be estimated via mass flux conservation at the surface, i.e.. 



Mace = 47r/7?^PpreVff 



(2) 



where the / is the filUng factor of the stellar photosphere 
heated by the accretion stream. The value of / was esti- 
mated to be 1.1% for the dipole magnetic geometry assumed 
above (see also Cranmer 2008). Assuming complete ioniza- 
tion (with 10% He by number), the preshock electron number 
density A'pre = 5.7 x 10^' cm~^. Because the shock is strong, 
we obtain an immediate postshock density A^post = 2.3 x 10^^ 
cm"^ (i.e., a factor of four increase across the shock, with the 
assumed adiabatic index of 5 /3). The postshock velocity Vpost 
is thus approximately 130 km s~' . 

As one proceeds deeper through the coohng zone, the den- 
sity increases and the velocity and temperature decrease, until 
a depth is reached at which there is no longer any influence 
from the shock. The "bottom" density phot can be estimated 
using ram pressure balance between the accretion stream and 
the unperturbed atmosphere. Using Equation 2 the pressure 
of the accretion stream is given by 



/9preV| ^ VffMacc 
2 87TfRl 



(3) 



(e.g., Hartmann et al. 1997; Calvet & GuUbring 1998). As- 
suming the accretion is stopped in the first few scale heights 
above the photosphere (at which radiative equilibrium drives 
the temperature to ~ (4/5)reff), the pressure of the atmosphere 
is 0.8/9botfe7eff . By equating the two pressures, the density can 
be found by solving for 



Pbot ' 



(4) 



with Teff = 4070 K for TW Hya (Batalha et al. 2002; see also 
Cranmer 2008). Converting to electron density, this corre- 
sponds to 2.6 X 10'^ cm"^, a factor of 1000 higher than the 
immediate postshock density. 
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Fig. 9. — Lower: Tg (dotted) and Ne (solid) as functions of depth in the 
postshock cooling column model, where corresponds to the position of the 
shock. The model is described in Section 4. Upper: Fractional contributions 
to the line intensities as labeled, as functions of depth in the model. 

Finally, we computed the physical depth and stracture of 

the postshock cooling zone using the analytic model of Feld- 
meier et al. (1997), which assumes that the radiative cooling 
rate Q^^a ~ NgA(T) has a power- law temperature dependence 
(A cx T"'/^). The thickness of the cooling zone is proportional 
to Vp^j,,/ ppost, and for the case of the TW Hya parameters used 
here, it has a value of 385 km. The spatial dependence of 
density, temperature, and velocity in the cooling zone is spec- 
ified by equations (8)-(ll) of Feldmeier et al. (1997), and 
the shapes of these curves resemble those of Giinther et al. 
(2007) — ^which were computed using a more detailed radia- 
tive cooling rate — reasonably well (see Fig. 9). 

Note that choosing the stronger mass accretion rate of 
BataUia et al. (2002) would have led to a truncation ra- 
dius closer to the star {r, = 2.2/?*), a lower freefall velocity 
(vff = 430 km s"'), and thus a lower postshock temperature 
(2.4 X 10* K). The densities would all have been higher than 
their counterparts above by about an order of magnitude, and 
the thickness of the cooling zone would have been smaller (32 
km). 

Using the temperature and density in the model, we com- 
pute the line intensities as functions of distance from the ac- 
cretion shock. Figure 9 shows line intensities through the re- 
gion for the strongest lines of Ne IX, O VIII, and O VII. We 
compare the integrated intensities to the observations. As- 
suming the abundances determined from Model D, and using 
Nh = 1.0 X 10^' cm"^, the predicted line intensities for the 
Ne IX, O VIII and O VII resonance lines are 0.25, 0.16, and 
0.40 times the observed fluxes, respectively. Higher abun- 
dances improve the agreement, as do lower column densities. 
Nevertheless, agreement within a factor of 4 seems quite rea- 
sonable, given the difficulty of obtaining absolute abundances 
and the complexity of the absorption. 

The predicted G-ratios give = 2.5 and 2.2 MK for Ne IX 
and O VII, respectively. The Ne calculated for Ne IX is 
2.6 X 10^^ cm~^, in excellent agreement with the observation; 
however, for O VII, the calculated A^^ is 4.1 x lO^'^ cm~^, seven 
times larger than observed. The O VII density is clearly dis- 
crepant with the standard model. 

5. ANALYSIS OF THE POSTSHOCK COOLING PLASMA 

The physical conditions predicted at the shock front are in 
excellent agreement with the observations; however, the pre- 
dictions for the postshock coohng plasma are in stark dis- 



TABLE 8 

Two-Region Model for Accretion Shock and 

Pns i'SiincK Cooling 



Model Parameters" 


Ne Te 


V Nh 


(cm-3) (MK) 


(cm') (cm"^) 


Region 1 6.0 X lO'^ 3.00 1 


5 X lO^** 1.0 X lO^i 


Region 2 2.0x10" 1.75 4 


5 X Hr" 4.0 X 10-" 


Predicted Unabsorbed Fluxes 


" from Region 1 


r i 


f Lya 


O 260. 166. 


3.6 512. 


Ne 447. 260. 


73.9 


Mg 3.6 1.3 


1.8 


Predicted Unabsorbed Fluxes 


" from Region 2 


r i 


f Lya 


O 124. 70.8 


36.7 25.1 


Ne 16.1 4.8 


10.1 


Mg 0.0 0.0 


0.0 


Ratio of Predicted Absorbed Fluxes" to Observed Fluxes" 


r i 


f Lya 


1.00 0.69 


0.53 0.69 


Ne 1.31 0.85= 


0.75 " 


Mg 1.28 1.14 


1.13 



" Abundances from Model D are 0.09, 1.04, and 0.18 relative 
to solar (Anders & Grevesse 1989), for oxygen, neon, and 
magnesium, respectively. 

* Fluxes at Earth are in units of 10~* ph cm~^ s~' . 
= Chen et al. (2006) is about 25% larger than APEC. 



agreement with the observations. Most problematic is that 
the electron density derived from O VII is lower than that for 
Ne IX, when the postshock gas should be compressing as it 
cools and slows and hence the density from O VII should be 
larger than from Ne IX. Furthermore, Nh determined from 
O VII is four times smaller than that determined from Ne IX, 
whereas in the standard postshock model O VII should be ob- 
served through the same or larger column density of interven- 
ing material. In this section, we explore the properties of the 
postshock cooling plasma based on our new results. 

Using the APEC code to calculate emissivities for a relevant 
grid of Te and Ne, we construct a set of two-component shock 
models for comparison with the diagnostic fines from the 
shock. The two components represent Region 1, the higher 
Te region near the shock front, and Region 2, the lower Te re- 
gion of postshock plasma. We do not discuss the hot coronal 
plasma here, but consider the diagnostic lines from ions that 
are formed below 3 MK, in particular O VII, O Vin, Ne IX, 
and Mg XI. Each component in the model has four parame- 
ters: Te, Ne, Nh, and volume V. The two regions experience 
different absorption with a factor of 2.5 larger Nh at higher 
Te, as required by the resonance series lines. The abundances 
are taken from the Model D fits. The predicted fluxes from 
the two regions are summed for comparison with the observa- 
tions. Table 8 gives the parameters for the best-fitting model, 
the predicted fluxes for each of the two regions, and the ra- 
tios of predicted to observed line fluxes. For O VIII, Ne IX, 
and Mg XI the predictions agree to better than 30% with the 
observed line fluxes. For O vn the predictions are low by a 
factor of less than 2, but the line ratios agree within 30%. The 
diagnostic line fluxes place tight constraints on the model. 

In this simple model, all of the emission from Mg XI and 
most of the emission from Ne IX come from Region 1, which 
we associate with the shock front. If we assume a cooUng 
length of 200 km, based on the accretion shock model de- 
scribed above (see Fig. 9), and use the derived emitting vol- 
ume of 1.5 x lO'^^ cm^, we obtain an area filUng factor for 
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Region 1 of 1.5% of the surface area of the star 

Region 2 has a larger volume than Region 1 by a factor of 
300 as determined primarily from O VII. If we assume that 
the cooling zone is the continuation of Region 1, we again 
obtain a length of about 200 km; however, this short length 
implies a filling factor larger than the stellar surface area and 
is not plausible. Instead, we use the scale height of 0.1 R^, 
for 1.75 MK plasma to estimate a filling factor of 6.8% for 
Region 2, a value in good agreement with estimates from the 
ultraviolet (Costa et al. 2000) and optical continuum veiling 
(Batalha et al. 2002). Scaling the mass as NeV implies that the 
mass of Region 2 is 30 times the mass of Region 1 . 

Romanova et al. (2004) have used three-dimensional mag- 
netohydrodynamic (MHD) simulations to show that the hot 
spots formed on the surface by accretion are inhomogeneous, 
with the lower density and temperature regions filling a larger 
surface area than the denser central region, seemingly sugges- 
tive of Region 2. However, the larger mass of Region 2 also 
implies a higher accretion rate than for Region 1 , and thus a 
larger density, in contradiction to our diagnostics. Therefore, 
we do not interpret Region 2 as a second accretion shock with 
a different accretion rate, but rather as a new accretion-driven 
phenomena. The impact from the accretion presumably pro- 
vides the energy to heat additional material in the stellar atmo- 
sphere to over 1 MK. Very recently, two-dimensional MHD 
simulations have shown that violent outflows of shock-heated 
material can propagate from the base of the accretion shock 
for cases with a high thermal to magnetic pressure ratio /3 
(Orlando et al. 2009). 

One of the key concepts of the standard one-dimensional 
accretion model suggests a simple picture (Fig. 10). The field 
line that truncates the disk separates the open magnetic field 
lines of the polar region from the closed magnetic field loops 
in the equatorial regions. The region heated by the shock (our 
Region 2) provides a large supply of ionized material avail- 
able to form the soft X-ray emitting corona presented here, a 
warm wind (Dupree et al. 2005), and soft X-ray jets (Giidel et 
al. 2008). This corona exists because magnetic loops confine 
it, but it is fed by the accretion process and is thus unlike stel- 
lar coronae in dwarf stars on the main sequence. Figure 10 
illustrates the shock and its surroundings. 

This model provides an explanation for the soft X-ray ex- 
cess discovered in the XMM-Newton Extended Survey of the 
Taurus Molecular Cloud (XEST; Gudel et al. 2007c; Telleschi 
et al. 2007). Grating spectra show that this excess emis- 
sion manifests itself in CTTS as enhanced O VII relative to 
O VIII as compared with weak-line T Tauri (WTTS) and main- 
sequence stars (Giidel & Telleschi 2007; Robrade & Schmitt 
2007). Our O VII measurements also demonstrate that the 
postshock cooling gas can, in some cases, dilute the signature 
of high Ne and explain why not all accreting systems show 
high densities. 

Having carefully distinguished between the accretion and 
corona in Section 3, we reconsider the assumption that these 
two components arise from different processes. Giidel & 
Telleschi (2007) find that the soft X-ray excess depends both 
on the level of magnetic activity as measured by the total (in- 
cluding coronal) Lx and on the presence of accretion. While 
we suggest that the role of accretion is to heat and ionize the 
surrounding stellar atmosphere to create a soft X-ray (post- 
shock cooling) plasma, it seems possible that this plasma 
could be further heated to 10 MK, either by the usual coro- 
nal heating (MHD) processes, or possibly by the accretion 
energy. Cranmer (2009) finds that in some systems, accre- 



Accretion-Fed Stellar Wind? 




Fig. 10. — Illustration of the shock and its surrounding environment. Re- 
gion 1 in the text corresponds to the shock and Region 2 to the postshock 
plasma. Figure courtesy of A. Szentgyorgyi. 

tion can provide the energy for coronal heating to such high 
temperatures. 

Studies of star forming regions have found that young 
stars show strong X-ray activity, with energetic flares and ex- 
tremely high X-ray emitting temperatures; however, the in- 
crease in activity with rotation rate breaks down for the ac- 
creting systems, leading to the conjecture that the stellar dy- 
namo might operate differently in these systems (Preibisch et 
al. 2005). We suggest that the accretion process may play a 
more prominent role than has generally been believed. While 
X-ray signatures of the accretion shock itself may not usually 
dominate the emission, the energy from accretion may in fact 
contribute to feeding and heating a new type of corona. 

6. SUMMARY AND CONCLUSIONS 

We have analyzed the deep Chandra HETG spectrum of 
TW Hya to assess the relative contributions of accretion and 
coronal emission. We summarize the most important of our 
findings: 

1 . The spectrum shows line and continuum emission from 
~ 10 MK coronal plasma, which may have a broad 
emission measure distribution. Unlike active main se- 
quence stars, the coronal EM distribution abruptly cuts 
off below ^ 5 MK, probably due to absorption by the 
accretion stream. 

2. The spectrum also shows emission lines from plasma 
at Tg = 2.5 MK, as predicted by standard models of the 
accretion shock. Tg and Ne diagnostics from He-like 
Ne IX are in excellent agreement with accretion models 
for TW Hya. 

3. We measure the flux of the O VII forbidden line in 
TW Hya for the first time to better than 3a, allowing a 
determination of electron density. This density is lower 
than Ne from Ne IX by a factor of four and lower than Ne 
from the standard postshock cooling model by a factor 
of seven, and leads to a new model of an accretion-fed 
corona. 
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4. The high Tg ions suggest a stellar corona. By con- 
trast, line ratio diagnostics from the lower Te ions, 
indicate two regions with different densities. In our 
model, the accretion shock itself has T^. = 3.0 MK and 
Ne = 6.0 X 10'^ cm"^, while the postshock cooling re- 
gion has Te = 1.75 MK and A^^ = 2.0 x 10" cm'^ The 
postshock plasma has 30 times more mass than the 
shock itself. The surface area filling factor of the shock 
is 1.5% while the postshock filling factor is approxi- 
mately 6.8%. 

5. Line ratio diagnostics require at least two different ab- 
sorbers, with higher column density (Nh) for higher 
charge state. Absorption by the accretion stream is 
larger by a factor of 2.5 for the shock than for the post- 
shock region. 

6. A range of models for the EM distribution and elemen- 
tal abundances provides acceptable fits to the global 
spectrum. While absorption, emission measure, and 
abundances are not independently determined by such 
methods, the high Ne/O abundance seems to be a robust 
result. 

7. From line profile analysis of Ne IX lines we determine 
a subsonic turbulent velocity of about 165 km s"'. 

This high resolution X-ray spectrum of TW Hya presents 
a rich set of diagnostic emission lines for characterizing the 
physical conditions of the high energy plasma and under- 
standing the dominant physical processes. The observations 



strongly support the model of an accretion shock producing 
a few MK gas at high density, as first suggested by Kastner 
et al. (2002) from an earlier HETG spectrum. The diagnos- 
tics also support the role of accretion in producing the soft X- 
ray excess at O VII previously discovered (Giidel & Telleschi 
2007; Robrade & Schmitt 2007). This excess emission re- 
quires that the accretion shock heat a large volume in the sur- 
rounding stellar atmosphere. In TW Hya the filling factor at 
O VII is larger than the filling factor of the shock at Ne IX, and 
is consistent with the filling factors of the ultraviolet and op- 
tical continua. Both open and closed magnetic field lines may 
emerge from this surrounding area to channel this ionized gas 
back into the corona. While it is not yet clear whether the 
energy from the impact of accretion at the star can heat a hot 
(10 MK) corona, accelerate a highly ionized stellar wind, or 
drive hot jets, we suggest that accretion-fed coronae may be 
ubiquitous in accreting young stars. 
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